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Abstract 

We discuss the propertíes of a large number N of one-dimensional (bounded) 
locally periòdic potential barriers in a finite interval. We show that the trans- 
mission coefficient, the scattering cross section a, and the resonances of a 
depend sensitively upon the ratio of the total spacing to the total barrier width. 
We also show that a time dependent wave packet passing through the system 
of potential barriers rapidly spreads and deforms, a criterion suggested by Za- 
slavsky for chaotic behaviour. Computing the spectrum by imposing (large) 
periòdic boundary conditions we find a Wigner type distributíon. We investi- 
gate also the S-matrix poles; many resonances occur for certain vàlues of the 
relative spacing between the barriers in the potential. 

PACS number(s): 03.65.Nk, 02.10.Yn, 05.45.Pq 
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I. INTRODUCTION 



Quantum systems with chaotic-like properties in the presence of tunneling through a 
single barrier have been studied recently [1]. The effect was attributed to the complexity 
of the wave function and its time dependence in the neighbourhood of the barrier, where 
Zaslavsky's criterion for the decrease of the Ehrenfest time [2] is expected to be satisfied. We 
discuss here several aspects of scattering from a one-dimensional locally periòdic potential 
barrier system [3,4] which is composed of a large number N of identical potential barriers 
densely arrayed along a finite section of the x axis. We discuss both cases of finite and 
infinite N. Parameters such as the transmission coefficient, the scattering cross section and 
the energy spectrum depend sensitively upon the ratio of the total spacing between the 
potential barriers to their total width. For example, it is shown in Section 2 for finite N 
and for both cases of e > v and v > e and in Sections 3, 4 and 5 for infinite N that as 
this ratio grows the transmission coefficient tends rapidly to unity. An approaching particle 
can, therefore, be transmitted unattenuated in its amplitude through these barriers without 
having to increase its energy, even when v » e (for the v > e case). Also, it has been 
shown, using the level statistics [5] of the energy spectrum of this dense system, that when 
this ratio increases the dense system appears to become chaotic-like in the sense of ref [1]. 
These chaotic-like characteristics emerge also, as will be seen in Section 4, when we study 
the passage of a Gaussian wave packet through the dense system. We show also that for 
both cases e > v and v > e the resonances of the scattering cross section depend strongly 
upon this ratio. 

Frishman and Gurvitz [6] have pointed out that the múltiple barrier structure is impor- 
tant to study. They analyse the finite múltiple well problem using a tight binding approxima- 
tion, and find a miniband structure (with some similarity to the Kronig-Penney spectrum) 
which may correspond to the property of ràpid approach to complete transmission that we 
find from certain vàlues of the ratio of spacing to barrier width denoted in the following by 
c. The exact solution that we study permits us to investigate the very sensitive dependence 
of the transmission phenomena, as well as the distortion of wave packets with high precision 
for both small and large number N of barriers restricted to a fixed total interval. We find a 
simple form for the limit N — > oo, which retains the very sensitive dependence on c. 

The configuration of a large density of one-dimensional barriers in a small range can 
occur in physical systems with planar defects such as successive evaporated layers or traps 
[7,8]; higher dimension analogs may also be realized. 

In section 2 we use a numerical model of this system where the number iV of potential 
barriers is finite. We discuss both cases of e > v and v > e. In Sections 3, 4 and 5 we discuss, 
using the transfer matrix method, the limit of an infinite number of potential barriers densely 
arrayed along a finite section of the x axis. The e > v case is discussed in Sections 3-4, and 
the v > e case in Section 5. In both cases we discuss the transmission probability, the 
scattering cross section, the poles of this cross section and the energy spectrum. We discuss 
also the level statistics of the energy spectrum of this dense system, and the properties of a 
time dependent wave packet that passes through it. 

We remark that although we study a system of N barriers that is locally periòdic, all of 
these barriers are contained, even for iV — > oo, in a finite interval. The system is therefore 
not equivalent to a Kronig-Penney type model [3] ; the emergence of the band like structure 
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as seen in Figure 3 appears to be due to the local (internal) periodicity but does not follow 
from global crystal translation symmetry, this periodicity is, moreover, not in the particle 
momentum, but in the total potential width. 

We believe that similar results may also be obtained for other similar systems like, for 
example, the N one dimensional periòdic potential wells in a finite interval. 



II. THE AN x AN MATRIX APPROACH FOR THE ENSEMBLE OF POTENTIAL 

BARRIERS 



The array of iV potential barriers discussed here is located along a finite section of the 
positive x axis beginning with the point x — 0. We denote the overall width of all the 
potential barriers by a, and the total width of all the interim spaces separating them by b. 
That is, in an array of N potential barriers the width of each one is and since in such an 
array there are (N — 1) separating spaces, the width of each one is rj^zï) • sketch of our 
array is shown in figure 1. To this array approaches from the negative half of the x axis a 
plane wave e lkx , where k = (-p 5 ) 5 . 

We consider both cases: e > v, and v > e, where v is the constant height of each potential 
barrier, and e is the energy of the coming wave function. We begin with the e > v case 
and write the following set of AN simultaneous linear equations obtained from the boundary 
conditions at the left and right hand sides of all the iV potential barriers [3,4] (see figure 1). 
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Re ik( a -^+b) + Se ^ k (^í+b) = Teiq c-^+b) + Ue - iq c-^+b) 
ikRe* 3 ^** - ikSe-^^^ = iqTe^^+V - iqUe***^** 

rp e iq(a+b) _|_ jj e -iq(a+b) _ ^ e ik(a+b) 

iqTe iq{a+b) - iqUe- iq{a+h) = ikZe ik{a+b \ 

where q = \J 2m ^í~ v ^ ~ . The former set can be written in a matrix form as Tx = (, where T 
is the square matrix with 4N rows and 4N columnns whose elements are given in the set 

(1). We denote by x the unknown vector with the 4N unknowns (A,B,C Z), and ( is 

the constant vector whose two first elements are —1 and —ik, and all its other elements are 
zero. As can be seen from the set (1) all the AN unknowns are obtained after dividing by 
the coefficient of the incoming wave so that \A\ 2 and \Z\ 2 (see the first two and the last two 
equations of the set (1)) are the coefncients of reflection and transmission respectively. In 
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order to calculate the value of \Z\ 2 we have to solve the 4iV simultaneous linear equations 
Tx = (. This can be done, especially for large N, only by numerical methods. We denote 
the constant length a + b of the x axis along which the system is arrayed by L, and define 
b = ac. Thus, we can express a and b in terms of L and c as follows 

T Lc 

b = -^~ (2) 



1 + c 1 + c 

The continuous curve in figure 2 shows the transmission coeficient \Z\ 2 as a function of 
c in the range 1 < c < 35 for N = 30, and the dashed curve is for N = 40. The other 
parameters are assigned the following vàlues: v = 100, e = 200, h — 1, m — |, L — 30. It 
is seen that the transmission coeficients, for both vàlues of N, tend to unity when c grows, 
but for the larger N, smaller vàlues of c suffice for the transmission coeficient to approach 
unity. That is, when the number of barríers íncreases, the approaching wave function passes 
unattenuated in its amplitude through these barriers even for relatively small vàlues of c. 
We note that we obtain the same result of a unity value for the transmission coeficient also 
when the restriction to a constant total length of the system is relaxed, as seen from Igure 
3. The dashed curve in Igure 3 shows the transmission coeficient as a function of the total 
width a when N = 60, and the continuous curve is for N = 120. For both curves we have 
assigned to b the value b = § . The potential v is 100 and the energy e is 200 as for Igure 
2. From the dashed graph, for N = 60, we see that the transmission coeficient has an 
almost constant periòdic pattern repeated as a function of a. In each one of these patterns 
the transmission coeficient oscillates near the value of 1, except at the beginning and end 
of each of these patterns where it drops to zero. In the continuous curve, for iV = 120, as 
in the dashed one, we have also a similar pattern repeated over the a axis, but this time 
the width of each such pattern is almost double. Checking the pattern of the dashed curve 
we see that the Irst drop of the transmission coeficient to zero occurs at a ps 9, whereas 
the corresponding drop in the continuous curve is at a ps 20. For N = 240 (this graph is 
not shown here) the transmission coeficient remains in the immediate neighbourhood of 1 
when 1 < a < 40. At about a = 40 this coeficient drops to 0, and then at a = 47 it rises 
to 1 and remains in the neighbourhood of 1 until a ps 87. Thus, as iV becomes larger the 
transmission coeficient appears to remain in the neighbourhood of unity for larger intervals 
of a, so that as N — > oo these intervals would appear to beco me inlnite in extent. In this 
limit one Inds agreement with the closed form we obtain in Eq (15). 

We, now, discuss the case v > e. In this case the set (1) has to be changed to take into 
account the tunneling intervals, that is, 
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(3) 



Ré k{ 
ikRe lk{ ~ 




In the set (3) k is the same as the k of the set (1) whereas q is q = ( ^ In (3) 
the wave functions inside the potential barriers contain real exponentials. This changes the 
former sinusoidal character of these wave functions (see the set (1)) to a hyperbolic one. 
The continuous curve in figure 4 shows the graph of the transmission coefficient \Z\ 2 as a 
function of c for N = 30, and the dashed curve is for N = 50. The potential v is 200, and 
the energy e = 180. The total length of the system and the range of c are the same as in 
figure 2, that is, L = 30, and 1 < c < 35. We see from both graphs that although v > e 
the transmission coefficient tends to unity as c grows, but this approach to unity is faster 
and for smaller vàlues of c, when N is larger. The same result is obtained if the condition 
of a constant total length of the system is relaxed as shown in figure 5, which shows the 
transmission coefficient as a function of the number N of potential barriers of the system. 
In this figure we take the total width a of all the potential barriers to be 8, and the total 
interval b to be | (c = |). The energy and the potential are assigned the vàlues of 200 and 
202 respectively. In this figure we see that the transmission coefficient has oscillating type 
behaviour when the number of barriers is small. At the larger vàlues of iV the transmission 
coefficient is in the close neighbourhood of unity. 

In summary, we see from this AN x AN matrix method applied to both cases of e > v 
and v > e, and for either a constant or variable total length of the system, that when the 
ratio c increases the transmission coefficient tends to 1, and when the number of potential 
barriers grows it tends to unity already at small vàlues of c. 

III. THE TRANSFER MATRIX METHOD FOR THE E > V CASE 

We discuss, now, the múltiple barrier system by the transfer matrix method [3,4,11], and 
in order to exploit its symmetry the dense array is assumed to be arranged between the 
points x = and x = where a and b has the same meaning as in the former section. 
We discuss first the e > v case. Using the terminology of Merzbacher [3] we can write the 
following transfer matrix equation which governs the behaviour of the bounded potential 
system. 



where A 2n +i and B<i n +\ are the amplitudes of the transmitted and reflected parts respectively 
of the wave function from the nth potential barrier. A is the coefficient of the initial wave 
that approaches the potential barrier system, and B is the coefficient of the reflected wave 
from the first barrier. is the product of three two dimensional matrices 
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The middle matrix T does not depend on n [12] and its components are given by 

ag .£ . ag .77 . ag 

Tu = cos( — ) + z- sm(— ), T 12 = 1- sm(-) 



(6) 
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As can be seen from Eq (5) the product of each neighbouring diagonal matrices M n M n _i 
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is constant for each n. That is, M n M n _x 
(4) as 



. Thus, we may write Equation 
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If we take the limit of a very large N, we obtain for the right hand side of the potential 
barrier system at the point x = ^ where n = N 



(9) 



We note that the last equation can be discussed from the eigenvalue point of view [3,4]. That 
is, by finding the appropriate eigenvalues from the suitable characteristic equation (a similar 
method has been applied to the finite N potential barrier system in [4]). In the following 
we adopt a more analytical and exact approach that yields the same results obtained from 
the former numerically-oriented method. 

The expression under the exponent N in Eq (9) can be written, using the set (6), in the 
limit of very large N, as 
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where 02 and 03 are the Standard Pauli matrices a 2 
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Using the 

relation lim fl _ ( . 00 (l + ^) n = e c , where c is some (possibly matrix-valued) constant we obtain 
from equations (9), (10) 
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We denote the two coeficients in the second exponent as 



making use of the relation 
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we can expand the second exponent on the right hand side of Eq (11) in a Taylor series. 
After collecting corresponding terms we obtain 
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Defining z — k(a + b), we , therefore, obtain 
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The corresponding expression for a single barrier of the same total width and location is [3] 
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One sees that the infernal structure of the múltiple barrier, in the limit of N — > 00, is 
different (they coincide only if b — 0). The determinant of the matrix on the right hand side 
of Eq (15) is unity. 

As we have seen, the equations (9)-(15) were concerned with expressing the amplitudes 
of the transmitted and reflected parts A 2 n+i, B 2 n+i of the wave function at the right hand 
side of the bounded system (at the point ^-^-) as functions of A , and -Bo at the left hand 
side of this system (at the point — ^ ^ ) . We can find these amplitudes at an arbitrary 
point — ^ < x < ^ by using the property of the system that it is an infinite sequence 
of potential barriers bounded at two sides, so the point x is associated with some barrier 
n and may be written as x = ±n(jj: + j^ï) = ±np, where p — + and n is in the 
range 1 < n < y. The potential cycle p can be expressed in terms of the total length L as 

Since x and L are finite numbers n must be infinite if 



so x 



P = f, so x = ±n w , orn-^ L 
N is. Thus, the amplitudes of the transmitted and reflected parts A 2n+1 , B 2n+ i of the wave 
function at the n-th potential barrier can be written as (compare with Eq (11)) 
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We have used Eqs (12) and the relation lim n ^ 00 (l + — ) 



e ClC2 , where c\ and C2 are 



arbitrary finite (possibly matrix-valued) constants. Now if we define j\ — di — ^ we 
obtain fl — d\ = = 4>\, where <p is given by Eq (13). Thus, we may use all the equations 
written before (for the right hand side of the dense system at the point x = ^) also for an 
arbitrary point — ^ < x < This result provides a closed form for the wave function in 
the potential region which we shall discuss further in a succeeding publication. 

(cos(ift) I i/ sinW ) 

Now, defining é K = , , we can find from Eq (15) the transmission proba- 

bility at the point x — - 2 
Thus, we obtain for this probability 
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where we have used Eqs (12). We see that the transmission probability reduces, when b — 0, 
to the known transmisson probability of the one potential barrier system which is located 
at the same place and exposed to the same wave function as the infinite potential barrier 
system [3,4], 
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We will see that the presence of a finite b in Eq (18) results in a new possibility for the 
transmisson probability to reach unity without having to increase the energy as in the one 
potential barrier system. Clearly, if a — * (no potential barrier), the transmission coefncient 
goes trivially to 1. Figure 6 shows a three dimensional graph of the transmission probability 
from Eq (18) as a function of the energy e and c. The total length is L = 70, v = 60, the 
range of e is 61 < e < 120, and that of c is 0.01 < c < 5. As expected the transmission 
probability tends to unity when the energy e grows, but as seen from the graph it tends 
faster, even critically, to the neighbourhood of unity (even for small e) as c increases through 
relatively small vàlues (i.e., this effect is not simply due to a — > 0), that is, as the total spacing 
becomes larger (see Eq (2)). This is in agreement with the results obtained in the previous 
section for the e > v case (see figures 2-3). 



IV. SCATTERING CROSS SECTION 

We now study the scattering cross section of the bounded potential barrier system. For 
this we use the S matrix which connects the outgoing waves A 2 n+i, and B at the two sides 
of the system to the ingoing ones -B2N+1 and A . That is, 
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Using the last equation together with equation (15) (we denote the two dimensional matrix 
from Eq (15) by Q where detQ = 1) one obtains for the four components of the matrix S 
[3] 
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Now, in order to find the phase shifts we have to find the eigenvalues A± from the following 
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According to the conventional phase shift theory [3,4] \± = e 2lS± , where S± are the phase 
shifts that correspond to the eigenvalues A±. From the last relations we obtain 
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The scattering amplitude is given by S — 1 = e 2í<5± — 1 = 2ie l5± sm(8±) = 2iriT±, and the 
cross section a± is then obtained as 

a± = 4tt 2 |T ± | 2 = 4sin 2 (5 ± ) = 2(1 - cos(25±)) (25) 

We see that a± — > for e — > oo, since d«O(^)->0 and ^ — > 1. It is further clear that for 
large e the period of oscillation grows; the oscillations with respect to e go as y/e, and the 
period, for which y/e + A = ^/e + 27T, is determined by A = 47T 2 + 4:7iy/e, therefore grows. 
Figure 7 shows a graph of function of the energy e (the same graph is obtained also 

for <7_). The total length L is 70, and a = 40, so b = 30. The potential v is taken to be 
70, and the range of e is 71 < e < 1000. One sees the increase in period on this graph, 
but the decrease in amplitude would not become visible until e » a 2 v 2 /4, i.e; for our case 
e » 2 • 10 6 . An interesting property of a± emerges when we relax the constraint of constant 
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L = a + b. In this case we find that the dependence of the cross section upon b is different 
from that upon a. That is, for the same value of a the cross sections <x±, as functions of 
the energy e, depend also on b only in a finite specific range that depends upon the value 
of a. For example, for v — 70 and a = 10, the cross section a + changes with b in the range 
of < b < 3.5a, and for b > 3.5a the change in a+ is so small that, as a function of e, can 
be considered constant. The same thing can be said also for <r_. We find that the cross 
sections depend on the total width a (for b fixed) in such a manner that the periods of a± 
are inversely proportional to a. That is, as the total width a of the potential barrier system 
grows the growth rate (with e) of the period of a± becomes smaller. 

We discuss now the energy level statistics [5] of the bounded dense array. To study 
this problem we use the S'-matrix and the boundary value conditions at two remote 
boundaries of the system. That is, using periòdic boundary conditions at the points 
\x\ = C, where C is much larger than the size L = a + b of the system, we obtain 
MN+i^ ikC = A e~ ikC , B 2 N+\e~ ikC = B e ikC . Thus, using the last two relations, and 
expressing the components of S in terms of those of Q (see Eq (21)), we write Eq (20) as 
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To obtain a non trivial solution for the vector 
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/ 2 sin 2 (0). , 2/sin(0) 



-) + 



(cos(2A;C + z) - cos(2^) cos(0)) - 2 sin(2A;C + z) cos(0)) = 



In ordcr to obtain the spectrum we solve, numerically, the last equation for the energies that 
satisfy both its real and imaginary parts. Obtaining these energies we use the unfolding 
procedure [5] to transform to the more appropriate energies from which we may obtain 
energy level distribution. We find that the distribution of the energy level spacings depends 
sensitively upon the value of c. For small c (small total interval b and large total width d) the 
relevant distribution is more of the Poisson distribution type [5] than of the chaotic Wigner 
one, whereas when c increases the corresponding distribution is more of the Wigner type 
than of the Poisson one. Figure 8 shows a histogram form of the level spacings distribution 
[5] of 102 energy levels obtained numerically for both cases of e > v and v > e. The potential 
height is here taken to be v — 120, and C = 90, L = 20 c = 19. The continuous curve is the 
chaotic Wigner distribution [5] as obtained from a random matrix model [9], and the dashed 
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one is the Poisson distribution [5]. One can see that the histogram-form curve resembles 
the chaotic Wigner one. The strong peak at 0.5 appears not consistent with the Poisson 
(dashed) curve. We remark that including levels only for e > v results in a distribution 
which is not clearly of Wigner type; if we select only v > e the distribution appears more 
clearly of Wigner type. We note that when we use Eq (27) for the v > e case, the /, d and 4> 
that must be substituted in this equation are not those defined by equations (12), (13), but 
those defined in the following section in Eq (34). 

In order to investigate further the properties of the bounded dense system we have passed 
a Gaussian wave packet through it and study its behaviour in the bounded potential region. 
We have used the complex packet 



where xq is the initial mean position of the packet in coordinate space, and po and wq are the 
initial mean momentum, and initial width (uncertainty) of the momentum respectively in p 
space. For our numerical simulations we have discretized space and time with a resolution 
of dx = | and dt — ^. This resolution ensures the condition dt < dx 2 which is necessary 
for a stable and steady performance of the numerical method used here [10]. For the other 
parameters we choose m = |, Wq = |, Xq = —10, c = 2.333 v — 2, and po — 3. The last 
two chosen vàlues ensure the condition of e > v. The dense system is arrayed between the 
points x = —10, and x = 10. The units we are using for length and time are therefore; 
x = ^jç- and t = (we take p to be momentum in units [mv] and w the dispersion in 
p). With this scale, we see that velocities in ^ are related to our parametric velocities 

^ cm = — ^ cm ■ During and after the passage of the wave packet through the potential 
region its initial Gaussian form is strongly deformed. The point-type curve in figure 9 shows 
the form of the density of the wave packet which evolves from Eq (28) when the number of 
potential barriers is N = 4, and the continuous graph is the form of this density for N = 150. 
The potential barriers arrayed between x = — 10 and x = +10 are not shown. Both curves 
in figure 9 are for the same time of t = 5.8. Comparing these curves we see that for large N 
the wave packet expands across the whole potential region more rapidly than the expansion 
for N = 4. It has been suggested by Zaslavsky [2] (see also discussion in Ref [1]) that this 
behaviour is characteristic of a system with classical Hamiltonian of chaotic type. Moreover, 
we see that the structure of the wave function for N = 150 is of a much higher degree of 
complexity. For larger vàlues of t (not shown) one sees that the forms of the transmitted 
and reflected waves are also of higher complexity. These results depend sensitively on the 
value of c. The effect is most pronounced when c is in the neighbourhood of 4, whereas, as 
c grows this effect diminishes until it completely disappears for very large c. Moreover, we 
obtain a significant transmission of the wave packet through the potential region even for 
the case of v > e, and the part of the wave packet that passes through this region increases 
as the number N or c or both of them grow. Note that we have found that the transmission 
of a plane wave grows rapidly with c, whereas it may be poorly transmitted when c is very 
small; in accordance with the wave packet behavior just described. 

We now study the problem of resonances associated with the dense system. In order to 
find them we find the resonances of the cross section tr± (see Eq (25)). Using Eqs (12), (13) 





(28) 
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and Eq (23) (or Eq (24)) we find that these are found at the vàlues of the energies e that 
satisfy the following equation 

1 + ^ = (29) 

Eliminating the total interval b (see Eq (2)) from the following equations and substituting 
in Eq (29) for d and from Eqs (12)-(13), and also for q and k we obtain 



sm 



{a + ac + cL) (e(fl + ac + cV) _ ^ + c))) = (3Q) 



2i i 

= ± — - ^ ye(a + ac + cL)(e(a + ac + cL) — va(l + c)) 

The last equation can be, of course, vàlid only if the energy e is complex. We denote this 
energy as e = e x + ie 2 - The components t\ and e 2 are found (see Appendix A) by solving 
the following two simultaneous equations. 

sin(r x 2 cos(^— tt— — )) cosh(r 1 2 sin(^— -^-^)) = ±r| cos(^-^-^^) (31) 

Zi Zi Zi 

cos(r 1 2 cos( ^ 1 — c ^ 7r ^ )) sinh(r 1 2 sin(^— ^^)) = ±r| sin( ^ 2 — c ^ ?r ^ ), (32) 
2 2 2 

where k — 0,1 and r±, cos(0i), sin(0!), r 2 , cos(02), and sin(0 2 ) are given respectively by 
Eqs (A3)-(A8) in Appendix A. We find numerically that there is no solution to equations 
(31), (32) for very large vàlues of c. The allowed range of c, for which these equations may 
be satisfied, depends upon the value of the total length of the system L; as L increases the 
allowed range of c expands. For all other vàlues of c outside these ranges we find no pole 
that satisfy the simultaneous equations (31), (32). We note that the poles are more freqüent 
at the middle sections of these ranges than at their ends. As noted the absolute vàlues of 
the complex energy must be greater than v since we deal here with the |e| > v case. It can 
be shown (see Appendix B) that as long as these absolute vàlues are not very much larger 
than the potential v the Eqs (31), (32) can be solved for a very large number of vàlues of ei 
and e 2 (dependent upon the vàlues of L and c). But when ei or both e\ and e 2 become very 
large these two equations have no solution for any value of L and c. 



V. THE V > E CASE 

We discuss, now, the v > e case. The matrix equations (4), (5) may also be used for the 
v > e case but the middle matrix T at the right hand side of Eq (5) has to be written as 

T n = cosh(^) + | sinh(^), T 12 = \ sinh(^) (33) 

Tj Q/Ç CLÇ ^ Q/Ç 

T 2 i = -- sinh(— ), T 22 = cosh(— ) - - sinh(— ) 
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k is y^, q is y/^j^, and £ and 77 are given by £ = -i V = f + f , r) = -i£ = & - f , 
where 77 and £ are from Eq (7). We can continue through the same steps as those of the 
e > v case and find that the equivalent of the /, d, and from Eqs (12), (13) are 

f = kb- ü -f, d=™f, 2 = (M - Úa 2 f = / 2 - d* (34) 
and the corresponding equations to (14)-(15) are 

eWa-sp**-^) = cos( V7^) + ^ - »^ 2 ) ^fp—fr) (35) 

V/ 2 - rf 2 







B2N+I 





e (cos0 + 



i/sin(<ft) ■ 



-tz dsin(<ft) 



— xe 

% . - 

é lz (cosó — i - sn ~^ ' 



A 
B 



(36) 



From the last two equations we can find the transmission probability for the v > e case (in 
an analogous way to the e > v case) 



*2iV+l ,2 



gÍ(K-z) 



^0 



(COS(0) 



_ ifsm(<j>) 



COS 2 (0) + 



/ 2 (sin 2 (0) 



1 + 



d 2 (sin 2 (0) 
^2 



(37) 



The e tK is the same as the e lK from the previous section (see the inline equation prior to Eq 
(18)) except that we substitute from Eq (34). Here, as for the e > v case, the last expression 
reduces, when b = 0, to the known transmission probability [3] for the one barrier located 
at the same place and exposed to the same wave function. 



A 



V cosh 2 (ag) + tíS^Ml 



(38) 



As for the e > v case the presence of a finite b yields a new possibility for the transmission 
probability to be 1 even when v » e. Figure 10 is a three dimensional surface of the 
transmission probability from Eq (37) as a function of the energy e and c. The range of e 
is 150 < e < 192, and that of c is 0.01 < c < 5 which is the same range as in figure 6. The 
potential v is 200. As for the e > v case the transmission probability tends to unity as c 
increases. 

We can continue, in a parallel way to the e > v case, and find the scattering cross section 
from the corresponding 5-matrix which is found from Eq (36) (We denote the matrix in this 

1 Q12 
Q21 1 

discussion for the e > v case) we can determine the corresponding eigenvalues À± which are 
the same as those found in Eq (22) except for the different /, d and 0. From these À± we 
can write equations corresponding to (23), (24), and from these equations the cross section 
à± may be obtained in an equivalent way to that of the e > v case. It is found that also 
here the period of à± becomes larger as e increases. As for the e > v case, if we release the 



equation by Q) to be S = -5^- 



From the Eq det(<S — XI) = (see the equivalent 
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condition of a constant L = a + b, we find that the dependence of the cross sections à±, 
as functions of the energy e, upon the total interval b is non-trivial only for a specific finite 
range which depends upon a. For example, for v = 140, e in the range 1 < e < 120, and 
a = 10 the cross section <r + , as a function of e, changes with b only in the range < b < 8.7a. 
For any other value of b > 8.7 a we obtain à + = 0. The range of dependence upon b becomes 
smaller as a grows. for example, when a grows from the former value of 10 to 15 the former 
range of b becomes smaller by almost a factor of 3, so that the new range in which cy u~ , as a 
function of e, changes with b is < b < 3.13a. That is, for any other value of b > 3.13a we 
obtain <r+ = 0. 

We note that although the periods of cr ± become larger as the energy e grows, the rate 
of growth is smaller compared to that of the e > v case, and as for the later case, the total 
width a and the rate of growth of the periods of à± are inversely proportional. That is, as 
a grows the growth rate of these periods becomes small. 

We discuss now the issue of resonances for the v > e case. We may use, for that matter, 
the equation corresponding to Eq (29) of the e > v case, except that we substitute the /, d 
and 4> from Eq (34) and also the q of the v > e case. We obtain the following equation 



sm 



a + ac + cL)_^ a + ^ + ^ _ ^ + ^ = ^ 



\ (1 + c) 2 

2i i 

= ± — — r\ eia + ac + cL)(e(a + ac + cL) — va(l + c)) 

av{± + c) v 

We now differentiate between two cases: The first is when v a(l + c) < e(a + ac + cL) in 
which case the last equation is identical to Eq (30), and so we can use the two simultaneous 
equations (31), (32) in order to find the real and imaginary parts e\ and €2 of the energy e. 
The essential difference between Eq (30) and Eq (39) (when va(l + c) < e(a + ac + cL)) is 
that in Eq (30) we have e > v, whereas here v > e. As for the e > v case, the allowed ranges 
of c depend upon the vàlues of L such that as L increases these ranges grow. 

We can prové, in an analogous manner to the e > v case (see Appendix B), that the 
energies that may be considered as poles of the cross sections cr ± can not assume very large 
vàlues (although here these energies have to satisfy the condition v > e). 

The second case is when we have in Eq (39) va(l + c) > eia + ac + cL), in which we 
obtain 



i sinh(, 



(a + ac + cL). . . r ... . Jn . 

' (av(l +c) - e{a + ac + cL))) = (40) 



\ (1 + c) 2 

2 / 

= =F — (Y+ j V e ( a + ac + cL)(va(l + c) — e(a + ac + cL)) 

This equation can be solved only for complex energies e = c\ + iei- As in the e > v case we 
eliminate the square roots and the complex character from both sides of Eq (40) by using the 
deMoivre theorem from Eq (Al) in Appendix A and the following hyperbolic sine addition 
formula [16] sinh(x + iy) = sinh(x) cos(y) + icosh(a;) sin(y). Thus, comparing separately the 
reals and imaginaries we obtain the following two simultaneous equations from which we try 
to find the real and imaginary parts of the energy 



14 



uf è Ai + ^nk . i . 0! + 27rÀ; i 2 + 2ttA; 
- cosn(7i cos( )) sm{r{ sm( )J = =Fr 2 cos( ) (41) 



sinh(rí cos( )) cos^ 2 sm( )) = ^r 2 2 sm( ), (42) 

V V V V 

where k — 0, 1, and f 1; cos(0i), sin^), f 2 , cos(0 2 ) and sin(0 2 ) are given respectively by 

fi = n, sin(0i) = -sin(0i), cos(0i) = -cos(0i) (43) 
f 2 = r 2 , sin(0 2 ) = cos(0 2 ), cos(0 2 ) = - sin(0 2 ) 

The variables n, cos(0i), sin(0i), r 2 , cos(0 2 ) and sin(0 2 ) are those of the e > v case and 
are given by equations (A3)-(A8) in Appendix A. It has been turned out, numerically, that 
there is no solutions to the two simultaneous equations (41), (42) that satisfy the condition 
of va(l + c) > e(a + ac + cL). 

We note that when we have released the condition of constant L we find (see the discus- 
sion on the cross-section à± before Eq (39)) that the cross sections à± become and remain 
zero for all vàlues of b that exceed some limiting value that depends upon the vàlues of a. 
We have, also, found that these limiting vàlues of b become smaller as a becomes larger. 
That is, in these cases the cross sections a± certainly have no poles. Now, when a becomes 
large the probability that the difference v a(l + c) — e(a + ac + cL) will be positive increases, 
and in this case, as we have just found, the cross sections à± have no poles in accordance 
with our discussion here. 

Summarizing the subject of poles in the last two sections we see that for the |e| > v case 
we find a very large number of different poles in large ranges of c, where the extent of these 
ranges depends upon the vàlues of L. When v > |e| the existence of poles depends upon the 
difference va(l + c) — \e\(a + ac + cl). That is, if this difference is negative then poles are 
found to the scattering cross sections à±, although in smaller ranges of c compared to the 
corresponding ranges of the |e| > v case. When the difference va(l + c) — |e|(a + ac + cL) 
is positive no pole of a± is found. 



VI. CONCLUDING REMARKS 



We have discussed the properties of a large number N of one-dimensional potential 
barriers arranged in a finite region of the x axis. We use both a AN x 4N matrix method for 
finite N and the transfer matrix method for the infinite array of these potential barriers along 
the finite region. We have discussed both cases of e > v and v > e and for both cases we 
found that the ratio of the total intervals between these potential barriers to their total width 
is an important parameter that determines the properties of the above mentioned variables. 
For example, when this ratio increases the transmission coefncient, for both cases of e > v 
and v > e, of the passing plane wave or wave packet tends to the unity value even when the 
initial energies of these waves are very small. A similar effect was found [8] in a classical 
diffusion system with a high density of imperfect traps for which the survival probability [7] 
of classical partides passing through it tends to unity when the interval between the traps 
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increases. Another system that was found [13,14] to demonstrate the same behaviour is the 
array of identical optical analyzers, such as Nicol prism [15], so that when the number of 
them, along a finite interval, becomes very large a beam of light passes through them with 
the same initial polarization and intensity it had before the passage. 

We have shown in this paper that a potential constructed of a large number of identical 
barriers can induce the type of behaviour as observed in the neighbourhood of tunneling 
barriers [1], interpreted in these references as chaotic-like. One may consider, as for the 
parallel drawn [1] between an unstable fixed point of the classical problem and quantum 
chaotic-like behaviour for the single barrier tunneling problem, a classical analog to the 
problem studied here. The repetitive potential in the bounded region, approached from 
above (for e > v), appears as an accumulation of unstable fixed points. The single wide 
barrier, on the other hand, is quasi-stable when approached from above; it is only the 
tunneling configuration in this case that has a strong analogy to the effect of a separatrix. 
For c large, when the potential barriers are relatively well separated, we see an apparent 
chaotic-like effect most strongly through the Wigner type level distribution, and when c is 
small the distribution moves toward Poisson type. 
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APPENDIX A 

Eqs (31)-(32) are obtained after eliminating the square roots and the complex nature 
from both sides of Eq (30). We do this by using the following two trigonometric relations 
[16]. 

/ / ,\ ■ ■ i ,wA i/ ,(p + 2irk. . . .(p + 2-Kk.. . . 

(r cos(0J + ir sm(0JJ« =r«(cos( J+ïsm( )), (Al) 

where n is any positive integer and k = 0, 1,2 . . .n — 1. 

sin(a ± ib) = sin(a) cos(iò) ± cos(a) sin(iò) = sin(a) cosh(ò) ± i cos(a) sinh(ò), (A2) 

Comparing separately the real and the imaginary parts of both sides we obtain the two 
simultaneous Eqs (31), (32) from which we can determine the components e\ and e 2 of the 
energies e that satisfy Eq (30). The six variables of Eqs (31), (32) that depends upon the 
cordinates (r, 0) are given by 



n 



(cl i (XC I CI J~j j ^ 

-((e^a + ac + cL) - va(l + c)) 2 + (e 2 {a + ac + cL)) 2 ) (A3) 



\ (i + cy 

cos(fc) = e^a + ac + cLÏ-avil + c) {M) 

J (ex(a + ac + cL) — av(l + c)) 2 + (e 2 (a + ac + cL)) 2 
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. / , x e 2 (g + ac + cL) 

sm(0i) = - (A5) 

J (ei(a + ac + cL) — av(l + c)) 2 + (e 2 (a + ac + cL)) 2 
r 2 = (A6) 



\Q((l -\- (xc + C.L) 2 

\ 4 4n n 4 — ((e 2 -ua(l + c) - 2eie 2 (a + ac + cL)) 2 + ((a + ac + cL)(ef - e\) - dva{\ + c)) 2 ) 

e 2 v a(l + c) — 2eie 2 (a + ac + cL) 



cos(0 2 ) 



e 2 f a(l + c) — 2eie 2 (a + ac + cL)) 2 + ((a + ac + cL)(e 2 — e 2 ) — eit>a(l + c)) 2 

(A7) 



. ,, . (a + ac + cL)(e? — eo) — eif a(l + c) 

sm(0 2 ) - - 



\J (e 2 t>a(l + c) — 2eie 2 (a + ac + cL)) 2 + ((a + ac + cL)(e\ — e 2 ) — e!fa(l + c)) 2 

(A8) 

APPENDIX B 

We show that there is no solution to Eqs (31)-(32) for very large vàlues of e x or of both 
ei and e 2 . In the first case we have t\ » e 2 and we obtain from equations (A3)-(A8) In 
Appendix A 

sm(<f>i) rí cos(0 2 ) 0, cos(0i) ss sin(0 2 ) rí 1, (Bl) 

ei(a + ac + cL) 2 4e 2 (a + ac + cL) 2 

n ~ (TTc] 2 ' r2 ~ a 2 v 2 (l + c) 2 

Using these approximations we can write the two simultaneous equations (31), (32) for the 
k = as (we note that the following two equations do not change their forms if this k 
assumes its second value of k — 1). 

. \ . NN . i /Tïki i 7T 7r/ci. , . 

sm(ri cos(7rfci)) cosh(ri sm(— — )) = ±r 2 cos(— H — — J (B2) 



cos(r 1 2 cos(7rfci)) sinh(r 1 2 sin(-^-)) = ±r 2 sin(— + 7r/ci) (B3) 

It can be seen that Eq (B3) is not satisfied for k\ = or an even k\. We, now, show that 
these equations are not satisfied for any uneven k\ either. For these k\ these two equations 
can be written as 

I I I 7T 

sin(r 1 2 ) cosh(ri ) = ±r| sin( — ) (B4) 



17 



I I I 7T 

cos(r 1 2 ) sinh(r 1 2 ) = ±r 2 cos(— ) (B5) 

Squaring the two sides of both equations we realize that their right sides are the same. So 
equating the left sides we obtain 

tan 2 (rf) = tanh 2 (rf ) (B6) 

i 

In order for the last equation to be vàlid the variable r{ must be small, but it is given that 
i 

ei is very large, so r{ must also be very large (see Eq (Bl)) and the equation (B6) can not 
be satisfied. Thus, when e\ » e 2 we find no poles of the cross section from (29). 

The same consequence is obtained also when both ei and e 2 are very large, so that we 
can write e x ~ e 2 . In this case we obtain 

sin(0i) cos(0i) ps -^=, cos(0 2 ) ~ —1, sin(0 2 ) ~ 0, 
V2 

and the corresponding two approximate simultaneous equations for this case are 

■ ( è r n , uf è • f n , nk \\ jl \ /7r + 27rA;i 

sin(r/ cos( — I — —)) cosn^r^ sm( — h — — )) = ±r 2 cos( ) (B7j 

8 4 8 4 2 

— 7T 7í/c^ i 7T 7ik\ — 7[k\ 

cos(r 1 2 cos(— H — — )) sinh(r 2 sin(— H — — )) = ±r 2 sin(— - ) (B8) 
8 4 8 4 2 

It can be seen that the first equation can not be solved for any ki (even or uneven). Thus, 
we see that very large energies can not be solutions of (31), (32). 
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FIGURES 



a/n a/r»b/(n-1) 3 aín t 2b/(n-1) (n-1)a/n»b a*b 

FIG. 1. The n potential barrier system. The approaching, transmitted and reflected waves are 
shown at right and left. 



Y The transmissíon coeffícient as a function of c=b/a. 
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FIG. 2. The continuous curve is the plot of the transmissíon coeficient for N = 30, and the 
dashed one is for N = 40, both as functions of c, and for the e > v case. Note that for the larger 
N the transmission coeffícient tends to unity for smaller vàlues of c. 



Y The transmission coelficients as functions of (he total width a 
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FIG. 3. The dashed curve is the plot of the transmission coeffícient for N = 60, and the 
continuous curve is for N = 120, both as functions of a, and for the e > v case. 
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Y The transmission coefficient as a function of c=b/a. 
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FIG. 4. The continuous curve is the plot of the transmission coefficient for N = 30, and the 
dashed one for N = 50, both as functions of c, and for the v > e case. Note that as in figure 2 the 
transmission coefficient for the larger N tends to unity at smaller vàlues of c. 



The transmission coefficient as a tunction ot the number N of potential barriers 




FIG. 5. The transmission coefficient as a function of N for the v > e case. The total width a 
is 8, the potential v is 202, and the energy e is 200. Note that this coefficient tends to unity as iV 
increases. 
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FIG. 6. A three dimensional graph of the transmission probability from Eq (18) for e > v as a 
function of c and the energy e. This graph is for a total system length of L = 70, v = 60, c in the 
range 0.01 < c < 5, and 61 < e < 120. Note that the transmission probability as a function of c 
tends to unity at a faster rate than as a function of e. 




FIG. 7. The graph of the cross section <7 + from Eq (25) as a function of the energy e. This 
graph is for a total length of L = 70, a = 40, v = 70, and for e in the range 71 < e < 1000. Note 
that the period of the graph becomes larger as the energy e grows. 
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FIG. 8. The histogram is the level spacing distribution as a function of s/D where s is the 
spacing between neighbouring levels and D is the mean spacing. The histogram is constructed 
from 102 energy levels found in the range 1 < e < 600, v is 120, C = 90, L = 20, and c = 19. The 
dashed curve is the Poisson distribution and the continuous graph is the Chaotic Wigner one. 
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FIG. 9. The continuous curve is the form of the wave packet for n = 150, and the point-type 
curve is the form for n = 4. Both curves are at time t = 5.8 in computer units. 
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FIG. 10. A three dimensional graph of the transmission probability from Eq (37) for v > e as 
a function of c and the energy e. This graph is for v = 200, L = 70, c in the range 0.01 < c < 5, 
and 150 < e < 195, We see from the figure that as the energy c grows above the value 0.2 the 
transmission probability jumps to 1. 
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